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Abstract 

The models for the internal structure of the newly found four-quark charmonium- 

like resonance Zc(3900) are discussed: the molecular model as well as the hadro- 

H ■ charmonium and tetraquark schemes. It is argued that it would be possible to resolve 

C^ ■ 

between these models by combining measurements of the quantum numbers of the 

resonance and of its decay rates into yet unseen channels TTf/;', vr/ic, prjc and into pairs 

of heavy mesons D*D and DD* . The models also predict different related four-quark 

states, which can be sought for in the existing and future data. 



The reported by BESIII pLj charged charmoniuni-hke peak Zc(3900) in the channel t^^J /ip 
is the newest addition to the 'collection' of known states related to heavy quarkonium that 
manifestly require the presence of two quarks and two antiquarks in their composition. 
Other hadrons of this type are the peaks in the i^^ifj' and vr^Xci spectra reported by Belle [21 
E] (although not confirmed by BABAR searches [H U]) and the bottomonium-like 'twin' 
resonances Zf,(10610) and Zfe(10650) produced in the decays T(5S') — t- ii'^Z^ and observed 
as peaks in the invariant mass spectra in the channels with bottomonium |6j (vr^T (nS") with 
n = 1,2,3 and 7i'^hf,{kP) with k = 1,2) and with pairs of B (B*) mesons [7J. Clearly, such 
states present a challenge for theoretical description of four-quark systems, and a better 
understanding of the internal workings of these and similar, yet unobserved, resonances may 
provide new insights into the strong dynamics of multiquark systems. 

The masses and the electric charge of the discussed resonances definitely require the 
quark composition to be Zq ~ QQud with Q standing for c or b, and the models discussed 
in the literature differ in the picture of clustering the quarks and antiquarks in this four- 
quark system that is used to somewhat organize and simplify the description. The so-far 
discussed models can be classified as follows. 

Hadronic molecules\^\. Heavy-light quark-antiquark pairs form heavy mesons, and the 
meson-antimeson pair moves at distances longer than the typical size of the meson. The 
mesons are interacting through exchange of light quarks and gluons, similar to nuclear force. 
Hadro- quarkonium^ [10]. The QQ pair forms a tightly bound system whose wave function 
is close to that of one of the heavy quarkonium states. The heavy quark pair is embedded in 
a spatially large excited state of light mesonic matter and interacts with it by a QCD analog 
of Van der Waals force. 

Tetraquarks^^. The pairs Qq and Qq form relatively tightly bound diquark and antidiquark, 
which interact by the gluonic color force. 

Clearly, other four-quark configurations are logically possible, e.g. a more uniform state 
where no significant pairing occurs. Furthermore, most likely all types of configuration con- 
sistent with the overall quantum numbers are, to some extent, present in the wave function 
and are quantum-mechanically mixed, and the difference between the discussed clustering 
models is in the assumed prevalent configuration with the other ones being considered as 
a relatively small admixture. Given this approximate classification, it is quite likely that 
the observed (and yet unobserved) four-quark states exhibit different type of the dominant 
behavior [9]. For instance, a description pjj of the ^^(10610) and Zfe(10650) resonances 
as being (dominantly) molecules made of respectively BB* {B*B) and B*B* pairs agrees 



with the relative strength and phase of the couphng of these particles to para- and ortho- 
bottomonium. On the other hand, the 'affinity', of some of the charmonium-like four-quark 
states to a particular state of charmonium (e.g. Z(4.43) — )■ ntp', Z\^2 -^ t^Xc\) indicates 
that they likely contain that particular state embedded in the light-matter excitation in the 
dominant part of the wave function. 

Currently it is not yet clear where the newly discovered resonance Zc(3900) fits in this 
classification, and an interpretation and further studies of this four-quark state can be quite 
instrumental in gaining understanding of multiquark heavy-light systems. The so far dis- 
cussed models describing this new peak include a DD* molecule or a cusp in the DD* 
spectrum [13] as well a molecular or tetraquark structure [H]. It has been concluded [13] 
that a molecular picture is likely preferred over the cusp hypothesis on the basis of the ob- 
served shape of the Zc(3900) peak, while in Ref. [H] detailed predictions of the tetraquark 
model for this resonance as well as expectation for related states are discussed and some 
of similar properties within the molecular model are also mentioned. In this paper I con- 
centrate on the behavior that should be expected in the molecular and hadro-charmonium 
models of the new state, and the ways of distinguishing by further experimental studies 
between the still open possibilities for its interpretation within the models. It will be argued 
that different models give distinctively different expected patterns of relative rates for the 
yet unobserved decays of the Zc(3900) resonance to the final states vr?/;', vr/ic, pric and DD* 
as well as different predictions for other related resonances. In what follows I first discuss 
the expected properties specific to the molecular model and then to the hadro-charmonium 
picture. 

In the molecular model the Zc(3900) is viewed as a resonance made of DD* and D*D 
pairs in the isovector state with positive G parity. (The I'^ = 1"*" assignment directly follows 
from the discovery mechanism for the resonance production: y(4260) — )■ 7rZc(3900).) If 
the heavy meson pair is in the S wave, as is also assumed in Refs. [131 113] , the spin-parity 
of the resonance is uniquely determined as J^ = 1^. In this case the pion in the decay 
y(4260) — )■ 7rZc(3900) is emitted in the 5* wavq^, and the chiral symmetry requires the 
amplitude of this process to be proportional to the pion energy E.,r (similarly to the behavior 



^The original paper [1] nientions a fit of tlie Zc resonance peak under tlie assumption tfiat tlie pion is 
emitted in tlie P-wave. It is however not clear whether there is an indication in the data that this process 
is a P-wave one, or that the assumption is ad hoc. Clearly, if the experiment would point toward a P-wave 
emission, the parity of the Zc resonance would have to be negative, and any discussion of it as an 5'-wave 
DD* molecule would be totally irrelevant. 



in the decays T{5S) -^ TcZb [T2]): 

A{Y -^ 7iZ,) (xE^{Y- Z) (1) 

with Y and Z being the polarization amphtudes for Y and Zc- In this picture the Zc(3900) 
resonance is a direct charmonium-hke analog of the bottonionium-like ^6(10610) resonance. 
Then a natural question arises of where is the analog of the higher Zb(10650) state? This 
leads to the expectation [H] that in this model there should be a similar 'twin' resonance 
Z'^ with the mass positioned relatively to the D*D* threshold similarly to positioning of the 
Zc(3900) relatively to the DD* threshold. The measured mass of the ^^(3900) is M{Zc) = 
(3899.0 ± 3.6 ± 4.9) MeV, so that the central value is 23.7 MeV above the D*+D^ threshold 
and 22.2 MeV above that for D^D*^. Placing the Z'^ resonance by the same amount above 
the ]J*^D*^ threshold gives its expected mass at approximately 4030 MeV. Assuming, as is 
the case for the ^^(lOGlO) and ^(,(10650) resonances, that the proportionality coefficients 
in the amplitudes, given by Eq.([T]), are approximately the same for the Z^ and Z^, one can 

r[y(4260) ^ vrz-] ^ r(z-^pzj/V.) ^ 

r[F(4260) ^ ttZJ ^ "-^^ ' ^^"^riz, ^ pij/^) ^ ^-^ ' ^^' 

so that the ratio of the combined transition rates through the two charmonium-like reso- 
nances can be estimated as 

i?r[y(4260) ^ 7r^Z,^(4030) ^ tt+tt- J/^] _ ^ ^^ 

5r[r(4260) -^ 7r±Z^T(3900) -^ tt+tt-J/^] ^ ' ' ^ > 

No peak of such significance is apparent in the data presented in Ref. pQ. However it is 
quite important that a dedicated experimental study of the presence of a peak near the mass 
4030 MeV similar to the Zc(3900) be done and an upper limit on its significance established. 
It is in principle possible that the ratio of the combined transition rates is somewhat smaller 
than the estimate in Eq.Q, e.g. due to a lager total width of the higher ^^(4030) resonance, 
which can be due to its coupling to the D*D channel. This coupling, which is suppressed 
by the heavy quark spin symmetry [12], can be enhanced for charmonium-like states in 
comparison with the behavior of the Z^ resonances due to lighter mass of the charmed 
quark. 

The S'-wave D*D molecular interpretation of the Zc(3900) resonance also implies distinc- 
tive properties of this state with regards to the total spin S of the cc quark pair. Namely, 
within this interpretation the spins of the heavy quark and antiquark are not correlated 
with each other, but rather with the corresponding light antiquark and quark. As a re- 
sult, in the molecular state the spin state of the cc pair is a mixture of S* = and S = 1. 



Specifically, in the I^{J^) = l^(l^) meson pair these two components are mixed with 
equal weight [12]. This behavior is well known for the Zi, resonances, which couple with ap- 
proximately equal strength to channels with ortho-bottomonium {TTT{nS)) and with para- 
bottomonium {7rhi,{kP)). Clearly, the same behavior should be expected of a molecular 
^2(3900), i.e. in addition to its decay into nJ/ip (and irip' discussed further in this text), it 
should have a comparable rate of decay into final states with para-charmonium: Zc — )■ vr/ic 
and Zc — ;• prjc- It should be noted that in this regard the molecular picture is somewhat 
similar to the tetraquark model, where the spin-correlated pairs are cq and cq, so that the 
state of the total spin of the cc pair is mixed. The ratio of the decay rates in the tetraquark 

model [11] is estimated as 

r(Z+ ^ p+r^J ^ 

which estimate does not look unreasonable in the molecular model as well. 

The known behavior of the Zi, resonances is that their pion transitions to excited T{2S) 
and T{3S) states are not suppressed (and rather enhanced) as compared to the transition to 
the lowest T(IS') bottomonium in spite of a significant kinematical enhancement of the latter 
one. This behavior is understood [15] in terms of larger overlap of the wave function of the 
bb quark pair in a spatially large molecule with spatially larger excited bottomonium states. 
The calculations [I5] based on modeling the wave functions of heavy quarkonium using 
the Cornell potential [16] are in a reasonable agreement with the data [7] on the relative 
rates of the pion transitions from Z^ resonances to various excitations of bottomonium. An 
application of the same approach to the pion transitions from Zc(3900) yields 

T[Zci390Q) ^ mP'] ^ 

r[z,(39oo) ^ ttJ/^] ^ ■ ^^ 

A similar, although a somewhat smaller estimate (about 0.3) for this ratio, is found in 
Ref. p^ within the tetraquark model. 

It can be argued that the interpretation of ^^(3900) as an S-wave D*D molecule possibly 
runs into difficulty related to its relatively high excitation energy, A ^ 23 MeV, over the 
threshold. Indeed, at such energy the characteristic momentum of the heavy mesons is 



p ~ y/MjjA ^ 200 MeV, where M^ is the mass of either of the mesons. Such momentum 
corresponds to a typical distance 1 fm between the mesons, which is uncomfortably close to 
the generally estimated range, where the mesons start to overlap, and can not be considered 
as individual particles. Also the interaction between the mesons should be quite contrived 
in order to explain a barrier that can hold an S wave resonance at about 23 MeV above the 
threshold. 



It is more natural if a resonance appears in a state with nonzero orbital momentum and is 
'held together' by the centrifugal barrier. The lowest nonzero orbital momentum corresponds 
to a P-wave motion of the heavy meson pair, which would imply negative parity for Zc(3900). 
In this case the options for J^ are 0~, 1^ and 2~, and the 1~ assumption can likely be 
discarded, since such state would have a large width due to decays into pairs of pseudoscalar 
D mesons, DD. The 0~ and 2~ cases are interesting in that in the I'^{J^) = 1"*"(0~) or 
I^{J^) = 1"*"(2~) state of D*D and DD* pairs the total spin of the cc pair is fixed, S* = 1. 
Indeed, the only combination of spin state of the heavy quark pair Sh and of the angular 
momentum of the rest (light) degrees of freedom J^ that has these quantum numbers is 
Iff (g) It for J^ = Q- and (1^ ® ij) © (1^ ® 2^) for J^ = 2". Thus if a future angular 
analysis finds the Zc(3900) resonance to be a J"^ = 0^ or a J'^ = 2~ state, one should 
expect in the molecular picture that the transitions from this resonance to states of para- 
charmonium, nhc and pr]c, are suppressed by the heavy quark spin symmetry relative to the 
transitions to ortho-charmonium, nJ/ip and irip'. The estimate in Eq.(|5]) for the relative 
strength of the latter two transitions should be applicable in this case as well. 

As is already mentioned, another possible interpretation of the newly found Zc{3900) 
is that it is dominantly a hadro-charmonium state, i.e. a tightly bound J/ip state of cc 
embedded in a light-quark excitated state. The observed decay Zc — t- nJ/ip is then a de- 
excitation of the light-quark matter. In this picture the Zc(3900) is tantalizingly similar to 
the resonance Z(4.43), which decays into irip', and has a very similar total width of about 
45 MeV. One can then view the latter resonance as a radial excitation of the cc pair over 
the Zc(3900) in essentially the same way as tp' is the radial excitation over the J/ip. The 
mass difference between Z(4.43) and ^^(3900) is approximately 535 MeV, which is by about 
55 MeV lower than the mass difference between ip' and J/ip, and one can speculate that this 
difference in the excitation energy can be attributed to the difference in the interaction with 
the light-quark 'environment' due to a larger spatial size of ip'. 

In the hadro-charmonium model the resonance Zc(3900) contains the cc pair in a pure 
5 = 1 state, so that the transitions to para-charmonium, Z^ — > vr/ic and Z^ — ?■ prj^ are 
expected to be suppressed. Furthermore, in as much as the cc pair has the wave function 
of J/ip (with possible slight distortions due to the interaction [9l [10]) the transition to 
tp', Zc — )■ 7cip', should be suppressed in comparison with the estimate in Eq.([5]). Another 
expected feature of the Zc(3900) resonance viewed as hadro-charmonium is that its decay 
into open charm, D*D and DD*, should be suppressed relative to the molecular case. In this 
respect the hadro-charmonium model is similar to the tetraquark scheme, where this rate is 



estimated [H] as T{Zc — ?• D*^D'^, D^D*^) ^ 4MeV, which accounts for only a small fraction 
of the total width of Zc(3900). In the molecular picture the significance of these decays can 
be gauged by the behavior of the Zh resonances, for which the 'dissociation' into heavy 
meson pairs constitutes (70-i-80)% of the total width [7], corresponding to the absolute rate 
of about 10 MeV or larger. One can expect that the decay into open flavor mesons should be 
enhanced for a molecular Zc(3900) due to its higher excitation energy above the threshold, 
so that these channels should account for a large, if not the major fraction of the total width. 

The simplest assumption about the quantum numbers of the Zc(3900) as hadro-charmonium 
is that it is the J ftp embedded in an S'-wave in a spinless excitation of the light-quark matter 
with the quantum numbers of a pion, i.e. J^ = 0~, so that the overall quantum numbers of 
the Zc(3900) are I^{J^) = 1+(1+). An assumption of a nonzero spin of the light-quark exci- 
tation and/or an orbital motion of the embedded J/tp would lead to a conclusion that there 
should also exist two or more states with nearby masses corresponding to a 'fine structure' 
due to the interaction of the spin of J/ip with the angular variables of the 'environment'. 
Since there appears to be no such structure in the data [1], it is reasonable to assume that 
the simplest arrangement of the hadro-charmonium embedding is realized in Zc- 

A distinctive prediction, stemming from a hadro-charmonium interpretation of Zc(3900), 
is that of an isovector four-quark resonance Wc, where the embedded J/tp is replaced with 
the rjc- In the limit of heavy quark spin symmetry the mass splitting between Zc and Wc 
should be the same as between J/ip and rjc, so that the expected mass of this lower resonance 
is 

M{Wc) ^ 3785 MeV , (6) 

and the expected dominant decay is Wc -^ tttjc with the same rate as Zc{3900) — )■ nJ/ip. It 
should be noted that Wc has to have the G-parity opposite to that of Zc, so that it cannot 
be produced in association with a pion in decays of y(4260). It can however be produced 
from higher 1 charmonium-like states in association with a p meson, e.g. e~^e~ — )■ p'^W^. 

Assuming the described simplest picture of the embedding for hadro-quarkonium, the 
quantum numbers of the Wc should be I'^{J^) = 1^(0+). A state with such quantum 
numbers and with mass given by Eq.(l6]) is certainly prone to a strong decay into DD pairs. 
However, in the hadro-charmonium picture the decay into open charm channels is expected 
to be inhibited [17] by the suppressed probability of the reconnection of the bindings between 
heavy and light quarks. It can therefore be expected that the total width of the Wc is not 
excessive as to prevent its observation in future experiments. 

It can be also noted that a resonance related to Zc{3900) with a lower mass is also 

6 



expected in the tetraquark scheme [H] and for a J^ = 1"^ molecule p!8]. In the former 
scheme the expected [Tl| mass "is about 100 MeV below " the Zf. resonance, which puts it 
distinctively higher in mass than the hadro-charmonium prediction ([6]). In the molecular 
scheme the J^ = O"*" molecule should be above the DD threshold by approximately the same 
amount as the J^ = 1~^ one is above the D*D threshold, which puts it at approximately 
3760 MeV, i.e. distinctively lower than given by Eq.([6]). It thus can be expected that a 
search for an I^ = 1~ charmonium-like resonance in the mass range 3750 -^ 3810 MeV will 
be helpful in resolving between the models of four-quark resonances. 

The main conclusion from the discussion presented in this paper is that a further ex- 
perimental study of the Zc(3900) resonance and related processes is vitally important for 
building an understanding of dynamics of multiquark heavy-light systems. It is argued that 
the most interesting aspects of such studies at the cm. energy corresponding to F (4260) are 

• establishing the spin and parity of ^^(3900); 

• a search for a peak around 4030 MeV in the nJ/ip invariant mass spectrum in the 
process F(4260) — > miJ/ip; 

• a measurement of the branching fraction for decays of ^^(3900) into heavy meson pairs, 

• a measurement of the rate of the decay Zc(3900) — )■ nip' relative to that of Zc(3900) — ?■ 

irJ/ip; 

• a search for the decays Zc{3900) — )■ Tihc and Zc — )■ prjc- 

Additionally at a higher energy of the e~^e~ beams a search for the hypothetical Wc resonance 
with the mass in the range 3750 -^ 3810 MeV can be performed using the process e'^e'^ — )■ 
pWc — )■ pTiric. As discussed above, quantitative data from these studies would allow to resolve 
between the models of the internal structure of the Zc(3900) resonance. 
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